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Abstract. Spherical accretion flows are simple enough for analytical study, by solution of the 
corresponding fluid dynamic equations. The solutions of stationary spherical flow are due to Bondi. 
The questions of the choice of a physical solution and of stability have been widely discussed. The 
answer to these questions is very dependent on the problem of boundary conditions, which vary 
according to whether the accretor is a compact object or a black hole. We introduce a particular, 
simple form of stationary spherical flow, namely, self-similar Bondi flow, as a case with physical 
interest in which analytic solutions for perturbations can be found. With suitable no matter-flux- 
perturbation boundary conditions, we will show that acoustic modes are stable in time and have no 
spatial instability at r = 0. Furthermore, their evolution eventually becomes ergodic-like and shows 
no trace of instability or of acquiring any remarkable pattern. 


RELATIVISTIC POTENTIAL FLUID FLOW 


We consider the adiabatic downfall of a perfect fluid onto a compact spherical body 
or a non-rotating black hole. In particular, we intend to analyse the stability of simple 
spherical, stationary flows (Bondi flows) [1]. General stability arguments have been 
given by Garlick [2] and Moncrief [3], but there are no exact solutions for perturbations, 
except in the WKB approximation. So we further restrict ourselves to self-similar flows, 
as a case amenable to analytic treatment and with physical interest. 

We begin with a summary of the theory of relativistic potential fluid flow and its linear 
preturbations, introducing the sound metric. Then we proceed to the Newtonian limit, 
sufficient for our purposes. In this limit, we obtain the self-similar Bondi flows and the 
perturbation equations. These equations can be solved in terms of Bessel functions. We 
study their initial and boundary problems, and so we draw conclusions on their stability. 

Let us consider the perfect fluid equations, namely, the energy-momentum tensor 
T^ v — (p + p)u^u v + pg ^ v , u^Ufx — —1. and thermodynamic equations h = (p+ 
p)/n , dp = n(dh — Tds). where n is the number density and h is the enthalpy per 
particle (we have = 0). The equations of motion are the conservation equations 
7^ v ;v = 0, ( nu = 0. Let us further consider isentropic solutions and potential flow 

[3], such that co^ v — {hu^)- v — ( hu v ) ;fl fulfills PaP^ft^v = 0, where P^ v — gp V + UpU v . 
Then we have ( 0 ^ v — 0 =>■ hu^ — \j/,p f° r some function yi. The equations of motion 
become 




where n is expressed in terms of h by the equation of state: n 
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Linear Perturbations and Sound Metric 


The linear perturbations of the equation for the scalar potential give the following 
scalar wave equation: 

v^w = o, 

where is the covariant derivative with respect to the sound metric 
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Therefore, causality in sound propagation is determined by (characteristics, etc), 
and the symmetries of are the ones common to g^ v and u - a . 

From the linear perturbation equation: 


V 


= 0, 


GfV _j_ 

» 2 


Sy^Sy K & VK - l - 8^ K<J 8y K 8y 0 


where ^ is the energy-momentum tensor of scalar waves. We consider stationary flow 
conserved energy: 

E = -2 I d 3 x(- detSf) 1 / 2 ^ = X -jd 3 x{- detSf) 1 / 2 [-& n (8y t ) 2 + & ij Sy/jSyj ]. 

In addition, we are interested in spherical flow =>■ conserved angular momentum. 


NEWTONIAN SELF-SIMILAR SPHERICAL FLOW 

When r>i?> GM/c 2 and \u l \ -C 1, potential flow boils down to v = Vt/r and [4] 
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This is an equation of non-homogeneous wave propagation. The law of conservation of 
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acoustic energy is simply ^ = V ■ W, where 
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and the acoustic energy flux current is 


W = ^d t Sy [\d t Sy - c 2 VSy + x ■ (\ - VSy)] = -d t Sy <5j, 


with j = pv the matter flux current. 



Basic flow: Bondi solutions 


For stationary spherical flow we have the radial mass conservation and Euler equa¬ 
tions 
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Integrating the first one, 4nr 2 pv = ,'/G\ (constant). We further use the polytropic equation 
of state p — JfopY. We define the accretion radius M = GM/(c s ) 2 „. and hence the non- 
dimensional variables 


r _ p _ v 
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Writing the equations of motion with these variables and solving for u : 

yj“V 2 + «(/ / "-l) =* _1 , 

where the adiabatic index n — l/(y — 1). 


Self-similar solutions 

Assume x«l^- ^-x~ 3 y~ 2 + ny l l n x — .r~ 3+2n z _2 " +nz— 1, where z — y l ! n x. If 

n = 3/2, we can solve for z(A) y — z(A) 3 / 2 x~ 3 l 2 , u = Az(A)~ 3 / 2 x _1//2 (power laws). 
Then, for n = 3/2 y= 5/3, 

p(r) = ar~ 3 / 2 , 

v(r) = /3 r -1 / 2 , 

c 2 (r) = ^ 2 ^a 2//3 r _1 = (J 2 r _1 , 

so the Mach number is given by ./# = v(r)/c(r) — /3/cr (constant). 


Linear Perturbations 


We try the separation of variables in spherical coordinates: 5yr(r, 0, <p,t) = 
P(r) Y[ m (6, (p ) So we obtain the radial differential equation 
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Note that it is not of Sturm-Liouville type. Its general solution is 

R(r) = r 1 / 4 g-'f"' 3/2 [Ci/ v (KT 3 / 2 ) +C 2 7- v (^ 3 / 2 )] , 
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We need appropriate boundary conditions at two radii. In the self-similar case, we 
must use r\ — 0 and r 2 ~ & —> °°, but we keep the latter finite to have a discrete spectrum. 
Since W = —d t Syf 5j, holding the mass flow, as done in Ref. [5], is equivalent to holding 
the energy flow. So, using the variable z — iff 3 / 2 , we impose 


z 1/6 (Ci[J v (z) + 6zJy(z)\ +C 2 [J-v{z) +6z/i v (z)]) =0 


at ri,r 2 , and so we obtain the K-spectrum. Remarkably, we have regularity at r\ — 0, 
namely, the physical quantities 8v r /v, Svq^/v, 8p/p stay finite. 


Evolution of radial eigen-functions 


Since our boundary problem is not of Sturm-Liouville type, we have no eigen¬ 
function orthogonality. Let us turn to the first order equations for radial perturbations: 


+-^d r [r 2 (Sp v +p Sv)] = 0, 

„2 


dt 


dSv 

dt 

-2; 


+ <9,- —5p+v<5v = 0 

p 


A = 


d(A-x) 
dr 

v(r) r 2 p(r) 
v(r) 


= iax, 


(rf 


r 2 p(r) 


o 1 


with x — {r 2 8 p,8v). If y —A-x, U = ^ j q J , then we get the orthogonality relation 
fy*-U-A~ l dr oc § nm . 

Then we express the initial condition y(r) as a sum of orthogonal modes: 


y{r) = Y, c nyn(r) = c 0 yo + 2Re ^c„y„(r) . 

n=—oo n= | 

Due to orthogonality, c n — {y n -,y)/ fn An) ■ The time evolution is given by c n (t) = 
c n e~ m,t . The norm (y, y) is invariant and, in fact, is proportional to the energy. The 
energy spectrum {|c„|~}“ =0 is invariant but the generic evolution of the correlation 
between the phases is to decrease with time, leading to quasi-ergo dicity and, therefore, 
(marginal) stability. 
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